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ABSTRACT 



The density fluctuation spectrum in the solar wind reveals a Kolmogorov-like 
scaling with a spectral slope of —5/3 in wavenumber space. The energy transfer pro- 
cess in the magnetized solar wind, characterized typically by MHD turbulence, over 
extended length-scales remains an unresolved paradox of modern turbulence theories, 
raising the question of how a compressible magnetofluid exhibits a turbulent spectrum 
that is characteristic of an incompressible hydrodynamic fluid. To address these ques- 
tions, we have undertaken three-dimensional time dependent numerical simulations of 
a compressible magnetohydrodynamic fluid describing super-Alfvenic, supersonic and 
strongly magnetized plasma fluid. It is shown that a Kolmogorov-like density spec- 
trum can develop by plasma motions that are dominated by Alfvenic cascades whereas 
compressive modes are dissipated. 

Key words: (magnctohydrodynamics) MHD, (Sun:) solar wind, Sun: magnetic fields, 
ISM: magnetic fields 



1 INTRODUCTION 

It is a curious observation that electron density fluctuations 
in the interstellar medium (ISM) ex hibit an omnidirectional 
Kolmogorov-like (|Kolmogorov 194~tl ) power spectrum k~ 5 ^ 3 
(or -11/3 spectral in dex in three dimensi o ns) over a 4 
to 5 decade range (jArmstrong et al 198ll : iHigdon 19841 : 
I Armstrong et al 19901 ). The solar wind plasma also pos- 
sesses density fluctuations that exhibit a Kolmogorov-like 
fc -5 / 3 spectrum (Goldstein et al 1995; Matthaeus & 
Brown 1988; Zank & Matthaeus 1990; Montgomery et 
al 1987; Lithwick & Goldreich 2001; Padoan & Nord- 
lund 1999; Shaikh & Zank 2007). Turbulent processes 
are believed to be re sponsible for the ob s erved density 
spectrum in the ISM IIArmstrong et al 198ll; IHigdon 1984 ; 
Higdon 198H; lArmstrong et al 1990l; lElmegreen 2004 ; 



iDickev et al 200ll ; 
eral MHD and 
simulations 
trum. Some 



iMinter fc Spangler 19961 ). Sev- 
hydrodynamic fluid turbulence 
show a Kolmogorov-like 5/3 spec- 
exam ple of which are ([Biskamp 2003; 



Shaikh & Zank 2007 



Lithwick fc Goldreich 2001 



Ghosh et al 19931 ; iKim fc Ryu 20051 ; iKritsuk 20071 ). An ex- 
haustive list of references describing a Kolmogorov-like 5/3 
spectrum in hydrodynamic fluid and magnetoplasma tur- 
bulence is how ever not possib l e to cite here . Some good re- 
views such as dMcComb 1990l; iLesieur 1990l; iBiskamp 20031 : 



Iar 



Matthaeus fc Brown 19881; 



iMatt: 
iMo 



Scalo fc Elmegre cn 2004; Elmcgrccn fc Scalo 2004)" and so- 

wind (|Goldstein et al 19951 ; 
Zank fc Matthaeus 199d; Montgomery et al. 19871 . 

Podesta et al 20071 ; IPodesta et al 20061 ; iPodesta et al 2008h . 
The Kolmogorov-like 5/3 spectrum is observed in many 
fluid, space and astrophysical plasmas as well. For in- 
stance, turbulent spectra in the ISM and galaxies are 
found t o exhibit a Kolmogorov-like scalin g in wavenumber 



IScalo fc Elmegreen 20041 ; lElmegreen fc Scalo 20041 ) discuss 
various possible scenarios of turbulence in general. While 
all of these works show the possiblity of a self-consistent 
energy exchange between widely disparate length-scales 
in the presence of waves, nonlinear structures, anisotropy, 
driving forces etc., the physical processes leading to a 
Kolmogorov-like turbulent density fluctuation spectrum 
is not yet fully understood. A great deal of attention has 
thus focused on understanding the evolution of MHD 
turbulence sp ect ra in the cont ext of the solar wind and ISM 
(IHigdon 19841; IHigdon 198rj ; iMatthaeus fc Brown 198~ 
Zank fc Matthaeus 19901 ; IZank fc Matthaeus 199. 



space (|Rov et al 200Sl; 
iRickett 20071; IWillett 2005: 



Haverkorn 20081; iRvu et al 2008t 



mcgrcen 



1999; 



Bavlv et al 19921; 



Montgomery et al. 198"! 



Lithwick fc Goldreich 200ll; IPadoan fc Nordlund 1999 : 
Shaikh fc Zank 20061 ; IShaikh fc Zank 20071 ). 
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Higdon (1984, 1986) interpret ed the observed density 
fluctuations l|Armstrong et al 198ll ) to be two-dimensional 
isobaric entropy variations in which temperature and density 
gradients are directed oppositely and both are orthogonal to 
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the local approximately unifo rm magnetic fie ld. Based on a 
pseudosound approximation (jLighthill 19521 ), Montgomery 
et al (1987) related density fluctuations and incompress- 
ible magnetohydrodynamics (MHD) velocity and magnetic 
field fluctuations. This approach, called a pseudosound ap- 
proximation, assumes that density fluctuations are propor- 
tional to the pressure fluctuations through the square of 
sound speed. The density perturbations in their model are 
therefore slaved to the incompressible magnetic field and 
the velocity fluctuations. This hypothesis was further con- 
trasted by Bayly et al. (1992) on the basis of their 2D com- 
pressible hydrodynamic simulations by demonstrating that 
a spectrum for density fluctuations can arise purely as a 
result of abondoning a barotropic equation of state with- 
out even requiring a magnetic field. Bayly's (1992) work 
ignores magnetic field effects from the outset and hence 
does not explain the influence and possible correlation of 
magnetic field and corresponding magnetized waves on the 
density fluctuation spectrum. The latter has been investi- 
gated in the slow solar wind plasma by Spangler & Spitler 
(2004) who suggest that there exists a strong correlation be- 
tween the density and magnetic field fluctuations. The pseu- 
dosound fluid description of compressibility, justifying the 
Montgomery et al. (1987) approach to the density-pressure 
relationship, was further extended by Matthaeus and Brown 
(1988) in the context of a compressible magnetofluid (MHD) 
plasma with a polytropic equation of state in the limit 
of a low plasma acoustic Mach number (Matthaeus and 
Brown, 1988). The theory, originally describing the gen- 
eration of acoustic density fluctuations by incompressible 
hydrodynamics (Lighthill, 1952), is based on a generaliza- 
tion of Klainerman and Majda's work (Klainerman and Ma- 
jda, 1981, 1982; Majda, 1984) and accounts for fluctua- 
tions associated with a low turbulent Mach number fluid, 
unlike purely incompressible MHD. Such a nontrivial fi- 
nite departure from the incompressibility state is termed 
a 'nearly incompressible' fluid description and is put for- 
ward to provide a possible explanation of the turbulent den- 
sity variations that are observed to exhibit a Kolmogorov- 
like power spectrum in the solar wind plasma (Montgomery 
et al. 1987; Matthaeus & Brown 1988, Zank & Matthaeus 
1990, 1993, Shaikh & Zank 2004a,b). In the context of the 
ISM, a comparative study of two-dimensional turbulence of 
self-gravitating supersonic MHD, hydrodynamic and Burg- 
ers turbulence by Scalo et al (1998) suggests a power law 
form of density fluctuations that is close to -1.7. Kim & 
Ryu (2005) report that the slope of the density power spec- 
tra in hydrodynamic turbulence becomes gradually shal- 
lower as the rms Mach number increases and it tends to- 
wards a Kolmogorov-like slope when the rms (or turbu- 
lent) Mach number is unity. The high resoultion simulations 
of Euler turbulence by Kritsuk et al (2007) suggest that 
the inertial range velocity scaling in the strongly compress- 
ible regime (with a spectral index close to -1.95) deviates 
substantially from the incompressible Kolmogorov 5/3-law. 
Kida & Orszag (1990) report that it is only the rotational 
component of the velocity field in driven hydrodynamic com- 
pressible fluid turbulence that exhibits spectra very close to 
that of the incompressible case even for a large Mach num- 
ber (close to unity). Lithwick & Goldreich (2001) show that 
density fluctuations in the ISM plasma are generated by 
entropy modes while Alfven waves lead to a Kolmogorov- 



like spectrum. This, ho wever, is not generic to all Alfveni c 
Mach numbers (Ma) jPassot fe Vazquez-Semadeni 20031 ) . 
While both slow and fast modes introduce density fluc- 
tuations at large Ma, only the slow mode dominates the 
dens ity-magnetic field anticorrelation at relatively small 
Ma l|Passot fc Vazquez-Semadeni 2 003). However, the slow 
mode is known to be strongly Landau damped in a colli- 
sionless plasma. The presence of a mean magnetic field in- 
troduces additional complexities in the energy cascade pro- 
cesses in MHD turbulence. For instance, the assumption of 
isotropy breaks down in the presence the mean or large scale 
magnetic field. Along the direction of this large scale mag- 
netic field Alfven waves suppress the parallel cascade (Irosh- 
nikov 1963; Kraichnan 1965). The resulting inertial range 
MHD spectrum is therefore thought to be flattened from 
£T 5/3 to £T 3/2 (Iroshnikov 1963; Kraichnan 1965). Within 
the paradigm of incompressible MHD turbulence, Goldreich 
& Sridhar (1995) proposed that parallel (fey) and perpen- 
dicular (fc±) modes are correlated by fen oc k 2 ^ 3 and energy 
in the perpendicular modes follows a Kolmogorov-like spec- 
trum when linear (along Bo) and nonlinear (across Bo) fre- 
quencies balance. By contrast, Shaikh & Zank (2007) argue 
that such balance is not obeyed identically by the entire in- 
ertial range modes, but only by a few modes that critically 
balance linear and nonlinear Alfvenic frequency. Boldyrev 
(2005) proposes a scale dependent anisotropic power law 
that appears to differ from the in-situ solar wind observa- 
tions (Podesta et al 2008). While in-situ spacecraft observa- 
tions of magnetic field fluctuations in solar wind plasma are 
shown to follow the Kolmogorov-like fc -5 ' 3 spectrum (Mont- 
gomery et al 1987; Zhou et al 1990; Goldstein et al 1995), 
the velocity field tends to show a close consistency with a 
fc~ 3/2 spectrum (Podesta et al 2006, 2007). The latter is 
contrasted by Roberts (2007) that the magnetic field and 
velocity in the solar wind do not evolve in the same way 
with helocentric distance. Based on Voyager observations, 
Roberts (2007a,b) argues that velocity spectrum relaxes to- 
wards a likely asymptotic state through spectral steepen- 
ing and acquires a spectral index of -5/3, finally mathch- 
ing the magnetic field spectrum. Roberts (2007b) further 
argues that -3/2 is accidental and transient, and that the 
-5/3 slope is the eventual state of all the fluctuations. In a 
comprehensive review, Bruno and Carbone (2005) and Vel- 
tri (1980) describe that low frequency solar wind velocity 
fluctuations closely follow a Kolmogorov-like spectrum, and 
the intermediate region (between high and low frequency) 
do not allow us to distinguish between a Kolmogorov spec- 
trum (-5/3) and a Kraichnan spectrum (-3/2). Bavassano et 
al (2005) describe that large scale velocity field fluctuations 
in the polar solar wind closely follow a Kolmogorov-like 5/3 
spectrum. 

The discrepancy in the magnetic and velocity field spec- 
tra continues to be an unresolved issue. Since our simulations 
tend to favor a Kolmogorov-like fc -5 / 3 velocity spectrum in 
which density spectrum (fc -5 / 3 ) closely follows the velocity 
spectrum through a passive convection, we support a fc -5 / 3 
Kolmogorov-like spectrum for the velocity field fluctuations 
in solar wind. What is notable in our work is the dissipa- 
tion of the high frequency component due to the damping 
(described in Section 2) of compressible plasma motion that 
suppresses the small scale and high frequency compressive 
turbulent modes. The MHD plasma therefore evolves to- 



The turbulent density spectrum in the solar wind plasma 3 



wards a nearly incompressible state and the density is con- 
vected passively by the velocity field to yield a k~ 5 ^ 3 spec- 
trum. 

Perhaps, the most striking point about the Kolmogorov- 
like fc -5 / 3 spectrum is its ubiquitous persistence in fluids 
and plasmas regardless of whether they are (in) compressible, 
(un)magnetized, (an)isotropic and (un)driven. The observed 
Kolmogorov-like density spectrum yields two paradoxes; (1) 
why does a compressible magnetized fluid behaves as though 
it were incompressible and umagnetized, and (2) Why do the 
density fluctuations, an apparently quintessential compres- 
sive characteristic of magnetized turbulence, yield a Kol- 
mogorov power law spectrum characteristic of incompress- 
ible hydrodynamic turbulence? These questions have to be 
answered if we are to address the origin of the 5/3 spec- 
trum in MHD turbulence in general and the ISM density 
power law spectrum, in particular. In this paper, we ad- 
dress these issues within the context of fully 3D simulations 
of compressible, anisotropic, driven Magnetohydrodynamics 
(MHD) turbulence to understand how and why a supersonic, 
super Alfvenic, and low plasma (3 (/3 the ratio of plasma 
pressure and magnetic pressure) compressible MHD fluid 
should exhibit a Kolmogorov-like wavenumber spectrum in 
density. We find a strong correlation between the intrin- 
sic MHD waves (i.e. AlfVen, fast & slow magnetoacoustic 
modes) and nonlinear inertial range turbulent cascades that 
suggests that nonlinear mode coupling interactions in com- 
pressible MHD turbulence tend to dissipate high frequency 
magnetoacoustic modes. Consquently, the inertial range cas- 
cade is governed predominantly by Alfvenic interactions that 
passively convect the density field to yield a Kolmogorov-like 
fc _5//3 spectrum. 

In section 2, we describe the governing equations, and 
critical assumptions of our 3D magnetohydrodynamic sim- 
ulations. Section 3 describes results from nonlinear fluid 
simulations of compressible, driven, anisotropic, homoge- 
neous, turbulent magnetohydrodynamic (MHD) plasma. 
Our simulations demonstrate that density, magnetic and ve- 
locity fields in MHD turbulence follow an omnidirectional 
Kolmogorov-like fc~ 5,/3 turbulent spectrum. The physical ar- 
guments explaining the evolution of a fc _5//3 spectrum are 
described in section 4 that deal primarily with turbulent 
damping of non-solenoidal velocity field fluctuations. In sec- 
tion 5, we outline our results for anisotropic cascades that 
result from the presence of a mean magnetic field in MHD 
turbulence. Finally, section 6 summarizes our major results. 



2 MHD MODEL 

Our underlying magnetohydrodynamic (MHD) model as- 
sumes that characteristic fluctuations in the magnetofluid 
plasma are initially isotropic, homogeneous, thermally equi- 
librated and turbulent. A large scale constant magnetic field 
is present and drives anisotropic turbulent cascades in an ini- 
tially isotropic spectral distribution in compressible MHD 
turbulence. The characteristic turbulent fluctuations in the 
plasma are assumed in our model to be much bigger than 
shocks or discontinuities. In other words, we ignore the in- 
fluence of shock on turbulent fluctuations. Our work incor- 
po rating the effect of shocks on turbu lent spectra is initiated 
in (|Zank et al 20071 ; IZank et al 2006r i. The boundary condi- 



tions are periodic, hence mode coupling interactions in the 
local spatial region are considered. 

The fluid model describing nonlinear turbulent pro- 
cesses in the magnetofluid plasma, in the presence of a back- 
ground magnetic field, can be cast into plasma density (p p ), 
velocity (U p ), magnetic field (B), pressure (P p ) components 
according to the conservative form 

^f + V-Q P = Q, (1) 
where, 
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1 rr2 Pp B 2 

Equations (1) are normalized by typical length £q and 
time io = £o/Va scales in our simulations such that 
V = l Q V,d/dt = t d/dt,tJ p = U p /Va,B = 
B/V A (^po) 1/2 ,P = P/poVlep = ep/p Vlp = p/p . 
The bars are removed from the normalized equations (1). 
Va = B /(^-Kp ) 1/2 is the Alfven speed 

The rhs in the momentum equation denotes a forcing 
functions (fivf(r, £)) that essentially influences the plasma 
momentum at the larger length scale in our simulation 
model. With the help of this function, we drive energy in 
the large scale eddies to sustain the magnetized turbulent 
interactions. In the absence of forcing, the turbulence con- 
tinues to decay freely. While the driving term modifies the 
momentum of plasma, we conserve density (since we ne- 
glect photoionization and recombination). The large-scale 
random driving of turbulence can correspond to external 
forces or instabilities for example fast and slow streams, 
merged interaction region etc in the solar wind, supernova 
explosions, stellar winds in the ISM, etc. The magnetic field 
evolution is governed by the usual induction equation, i.e. 
Eq. {l}, and obeys the frozen-in-field theorem unless dis- 
sipative mechanism introduces small-scale damping. Note 
carefully that MHD plasma momentrum equation contains 
dissipative terms on the rhs. It is the term in the momen- 
tum equation (i.e. d(p p \J p )/dt ■ ■ ■) that is proportional to 
^iV 2 Up + r;V(V'Up). The latter (along with the other terms 
in the equation) is divided by the density p p field during the 
evolution to calculate the velocity field. 



3 NONLINEAR SIMULATIONS 

We have developed a fully three dimensional compressible 
MHD code to study the nonlinear mode coupling interac- 
tion in the context of compressible MHD turbulence. De- 
tails of our code ar e described in (|Shaikh fc Zank 20061 ; 
IShaikh fc Zank 20071 ). In the simulations, all the fluctua- 
tions are initialized isotropically with random phases and 
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amplitudes in Fourier space and an initial spectral shape 
close to k~ 2 (k is the Fourier mode, which is normalized 
to the characteristic turbulent length-scale Zo). We have 
carried out the simulations for both decaying and driven- 
dissipative cases with and without external magnetic field 
i?o- Since we are interested in a local region of the solar 
wind magnetofluid plasma, the computational domain em- 
ploys a normalized three-dimensional periodic box of vol- 
ume 7r . Other parameters are 7 — 5/3, j3 = 0.1 — 2.0, Ma = 



1.0 - 2.0, M a = 1.0 - 2.0, 77 = n = 10" 



10" 



,k f < 15.0, 



where 7, Ma, M 3 ,r], n and kf are respectively the ratio of 
specific heats, Alfven Mach number, sonic Mach number, 
viscosity, magnetic diffusion and energy injection modes. 
Our MHD model does not include (photo)ionization and ra- 
diation terms. Hence source or sink terms corresponding to 
the density fluctuations are not included in our MHD model. 
Additionally, the localized dissipation is effective in our sim- 
ulations at the small-scales where it damps the plasma mo- 
tions. The small scale dissipation in the local interstellar 
medium or solar wind may r esult from e.g. radiative cool- 
ing or ion-neutral collisions (Spa ngler 19911) . Accordingly, 
the small-scale dissipation in our model corresponds to col- 
lisional or viscous effects and is associated with the small- 
scale damping that is responsible only for the cascade of 
large-scale energy into the smaller scales thereby producing 
a well-defined inertial range turbulent spectrum. By con- 
trast, the large-scales and the inertial range turbulent fluctu- 
ations remain unaffected by direct dissipation of the smaller 
scales. 

The initial kinetic and magnetic energies are equi- 
partitioned between the velocity and the magnetic fields. 
The latter helps treat the transverse or shear Alfven 
and fast/slow magnetosonic waves on an equal footing, at 
least during the early phase of the simulations. Magne- 
toplasma turbulence evolves under the action of nonlin- 
ear interactions in that larger eddies transfer their energy 
to smaller ones th r ough a forward cascade. Accordi ng to 
jKolmogorov 194ll ; llroshnikov 19631 ; iKraichnan 19651 ), the 
cascade of spectral energy is mediated by a local interac- 
tion amongst the neighboring Fourier modes that contin- 
ues until the energy in the smallest turbulent eddies is dis- 
sipated due to the finite Reynolds number. This leads to 
the damping of small scale motions as well. This results in 
a net decay of the turbulent sonic Mach number M s as- 
sociated with the large scale fluctuations. If turbulence is 
not driven at large scales, the turbulent sonic Mach num- 
ber continues to decay fro m a supersonic (M„ > 1) to a 
subsonic (M s < 1) regime ([Shaikh fc Zank 20061 ). This in- 
dicates that nonlinear cascades predominantly cause super- 
sonic MHD plasma fluctuations to become subsonic. In our 
decaying turbulence simulation, the large-scale energy sim- 
ply migrates towards the smaller scales by virtue of nonlin- 
ear cascades in the inertial range and is dissipated at the 
smallest turbulent length-scales. On the other hand, spec- 
tral transfer in driven turbulence follows a similar cascade 
process as in the decaying turbulence case. However, the 
inertial range spectrum in the latter is maintained by a 
large scale forcing at k < 5. The spectral transfer of tur- 
bulent energy in the neighboring Fourier modes in glob- 
ally isotropic and homogeneous hydrodynamic and magne- 
tohydrodynamic tu rbulence is the widely accepted paradigm 
|Kolmogorov 194lh that leads to Kolmogorov-like energy 



Figure 1. (Left) Velocity fluctuations are dominated by shear 
Alfvenic motion and thus exhibit a Kolmogorov-like k~ 5 / 3 spec- 
trum. The middle curve shows the magnetic field spectrum. Den- 
sity fluctuations are passively convected by the nearly incompress- 
ible shear Alfvenic motion and follow a similar spectrum in the 
inertial range. The numerical resolution in 3D is 512 3 . (Right) 
The evolution of Alfvenic (fc^) and fast/slow magnetosonic (&Afs) 
modes demonstrates that the spectral cascades are dominated by 
Alfvenic modes. 



spectra, while (Iroshnik ov 19631 ; IKraichnan 19651 ) describe 
turbulent spectra in the presence of a mean or local Bo. 
The most striking effect, however, to emerge from the de- 
cay of the turbulent sonic Mach number is that the den- 
sity fluctuations begin to scale quadratically with the sub- 
sonic turbulent Mach number as soon as the compressive 
plasma enters the subsonic regime, i.e. 8p ~ C(M S 2 ) when 
M s < 1. This was demonstrate d in our 3D simulat ions 
of a compressible MHD plasma l|Shaikh fc Zank 20061 ). It 
signifies an essentially weak compressibility in the magne- 
toplas ma, and is consistent with a nearly incompressible 
state (uMatthaeus fc Brown 198sl; IZank fc Matthaeus 199fj ; 



IZank fc Matthaeus 19931 ; IShaikhfc Zank 20061 ). 

In the context of the magnetoplasma being nearly in- 
compressible, the density fluctations exhibit a weak com- 
pressibility in the gas and are convected predominantly pas- 
sively in the background incompressible fluid flow field. This 
hypothesis can be verified straightforwardly by investigat- 
ing the density spectrum which should be slaved to the in- 
compressible velocity spectrum. This is shown in Fig. © 
which illustrates that the density fluctuations follow the ve- 
locity fluctuations in the inertial regime over the long time 
(several Alfven transit time) evolution of MHD turbulence. 
The evolution of compressible magnetoplasma from a(n) 
(initial) supersonic to a subsonic or nearly incompressible 
regime is gradual and it results in the density field follow- 
ing the velocity fluctuations. In the subsonic regime, com- 
pressibility weakens substantially so that density fluctua- 
tions are advected only passively. A passively convected fluid 
exhibits a similar inertial range sp ectra as that of its back- 
ground flow field (|McComb 19901 ). Likewise, subsonic den- 
sity fluctuations in our simulations exhibit a Kolmogorov- 
like fc -5 / 3 spectrum similar to the background velocity flue- 
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tuations in the inertial range. This, we believe, provides a 
plausible explanation for the Kolmogorov-like density spec- 
trum observed in MHD turbulence i.e. they are convected 
passively in a field of nearly incompressible velocity fluc- 
tuations and acquire identical spectral features [as shown 
in Fig. fT}]. The passive scalar evolution of the density 
fluctuations is associated essentially with incompressiblity 
and can be understood directly from the continuity equa- 
tion as follows. Expressing the fluid continuity equation as 
(<9t+U-V) hip = —V-U, where the rhs represents compress- 
iblity of the velocity fluctuations, shows that the density field 
is advected passively when the velocity field of the fluid is 
nearly incompressible with V-U ~ 0. The theoretical ba- 
sis illustrating the nonlinear damping of the non-solenoidal 
component of velocity field, i.e. V • U — > 0, is described 
quantitatively in the following section. 



4 TURBULENT EVOLUTION OF 

NON-SOLENOIDAL VELOCITY FIELD 

Understanding, how an initially non-solenoidal velocity field 
evolves towards a solenoidal field is important as it explains 
the evolution of a compressible MHD magnetoplasma from a 
supersonic to a subsonic or nearly incompressible state that 
yields a passively advected Kolmogorov-like k~ 5 ^ 3 density 
spectrum. Any given velocity field can be decomposed into 
solenoidal and non-solenoidal components or, equivalently, 
longitudinal and transverse components, which we shall here 
refer to as fast/slow magnetosonic and Alfven components, 
respectively. The simulations show that the amplitudes of 
fast/slow magnetosonic components described by the quan- 
tity /cms decay more rapidly than the amplitudes of the 
Alfvenic components. Consequently, they will be inefficient 
in cascading the corresponding inertial range spectral en- 
ergy. Hence nonlinear interactions, in the inertial range, are 
governed predominantly by non-dissipative Alfvenic modes 
(/ca) that survive collisional damping in compressible MHD 
turbulence. This is quantitatively demonstrated in Fig. 
{TJ [right panel]. 

The damping of a non-solenoidal velocity component, 
in part, explains the observed Kolmogorov-like k~ 5 ^ 3 in our 
simulations. To this end, it is essential to distinguish the 
Alfvenic and non-Alfvenic, i.e. corresponding to the com- 
pressional or due to slow and fast magnetosonic modes, con- 
tributions to the turbulent velocity fluctuations. To identify 
the distinctive role of Alfvenic and fast/slow (or compres- 
sional) MHD modes, we introduce diagnostics that distin- 
guish the energies corresponding to Alfvenic and slow/fast 
magnetosonic modes. Since the Alfvenic fluctuations are 
transverse, the propagation wave vector is orthogonal to the 
velocity field fluctuations i.e. k _L U, and the average spec- 
tral energy contained in these (shear Alfvenic modes) fluc- 
tuations can be computed as 



<Mt)> 



'E k l*xu k | : 



£ k |u k 



The above relationship leads to a finite spectral contribution 
from the |k x Uk| characteristic turbulent Alfvenic modes. 
On the other hand, fast /slow magnetosonic modes propagate 
longitudinally along the velocity field fluctuations, i.e. k || U 



Figure 2. Turbulent energy associated with the characteristic 
inertial range modes corresponding to |k X Uk| and k || U fluctu- 
ations is shown. Our simulations show that the velocity fluctua- 
tions are dominated by shear Alfvenic motion whose contribution, 
corresponding to the curve represented by |J7(fc)c ur iu| 2 , is more 
than an order larger than that corresponding to the fast/slow 
magnetosonic modes and is shown by | C/ ( A;) oivU | 2 - This result is 
consistent with Fig 1 (right panel) that depicts the mode coupling 
evolution of the two MHD modes. 



and thus carry a finite component of energy corresponding 
only to the ik • Uk part of the velocity field, which can be 
determined from the following relationship 



(k M s(t)) 



'E k l*-Ukl ; 



The expression of /cms essentially describes the modal en- 
ergy contained in the non-solenoidal component of the MHD 
turbulent modes. 

The quantative evolution of the characteristic modes 
corresponding to the Alfvenic /ca and slow/fast compres- 
sional magnetosonic /cms modes is depicted in Fig. 1 (right 
panel). Although the modal energies in /ca and /cms modes 
are identical initially, a disparity in the cascade rate de- 
velops, and the energy in longitudinal (or compressional) 
fluctuations associated with the non-solenoidal velocity field 
decays far more rapidly than the energy in the Alfvenic 
modes. The Alfvenic modes, after a modest initial decay, 
sustain the energy cascade processes by actively transfer- 
ring spectral power amongst various Fourier modes. By 
contrast, the fast/slow magnetosonic modes (kMs) progres- 
sively weaken and suppress their corresponding spectral con- 
tribution in the turbulent energy cascades. The difference 
in the cascades corresponding to /ca and /cms modes per- 
sists even at long times. The /cms mode represents collec- 
tively a dynamical evolution of small-scale fast plus slow 
magnetosonic cascades. The physical implication, however, 
that emerges from Fig. 1 is that the fast/slow magnetosonic 
modes (/cms) do not contribute efficiently to the spectral 
transfer process, and that the cascades are governed pre- 
dominantly by non-dissipative Alfvenic modes that survive 
the collisional damping in compressible MHD turbulence. 
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Correspondingly, the turbulent energy associated with the 
Alfvenic modes makes a dominant contribution to the ve- 
locity fluctuation spectrum when compared to the magne- 
tosonic modes. We have made measurements of the tur- 
bulent energy in the Alfvenic and fast/slow magnetosonic 
modes to quantify their respective contributions to the ve- 
locity field fluctuation spectrum. This is shown in Fig. 
(2). In Fig (2), |[/(fc) Cur iu| 2 = EJ k x U klVE fc fc2 and 
\U(k) DivV \ 2 = EJ k • U k| 2 /E fc fc2 - Clearl y- the ener gy 
contribution by Alfvenic modes (parallel to Bo and or- 
thogonal to the velocity field) is more than 10 times that 
of the fast /slow magnetosonic modes (parallel to the ve- 
locity field). This clarifies that it is the predominance of 
Alfvenic modes (Fig (1) & (2)) in inertial range cascades 
that primarily lead to a Kolmogorov-like spectrum. Damp- 
ing of Alfven wa ves is possible by ion-neutral collisions as 
pointed out by dKulsrud fe Pearce 19691 ; iBalsara 19961 ) in 
the context of molec ular clouds [for more references, see 
( Shai kh" &c Zank 20081 )]. In our present simulations, we do 
not include ion-neutral damping as we focus mainly on the 
solar wind plasma. This nonetheless suggests that because 
of the decay of the fast /slow magnetosonic modes in com- 
pressible MHD plasmas, supersonic turbulent motions be- 
come dominated by subsonic motions and the nonlinear in- 
teractions are sustained primarily by Alfvenic modes there- 
after; the latter being incompressible. One of the implica- 
tions of the turbulent damping of non-solenoidal velocity 
field (ik • Uk) in an MHD fluid is that compressible modes 
make negligible or no contribution to the inertial range en- 
ergy cascade. The cascade is thus determined primarily by 
the incompressible Alfven modes that passively convect the 
density fluctuations. This point is further consistent with 
the MHD fluid continuity equation which in k space reads 
as follows. 



*) 



-ik ■ U k . 



^+i^(5(k + k')U k -k'J in P (k, 

The nonlinear mode coupling interations associated with 
the Dirac delta function 8 are finite only for those inter- 
actions that obey the Fourier diad k + k' = in spec- 
tral space. It follows from our simulations that the tur- 
bulent damping of the non-solenoidal velocity field on the 
rhs of the continuity equation makes an insignificant contri- 
bution to the inertial range energy cascade. The nonlinear 
mode coupling interactions in MHD turbulence are there- 
fore dominated by convective transport that leads to a pas- 
sive convection of density fluctuations. The density fluctu- 
ations subsequently follow the inertial range spectrum and 
are identical to the background nearly incompressible veloc- 
ity field and thus have a Kolmogorov-like fc~ 5 ^ 3 spectrum 
(|Kolmogorov 194ll ; iMcComb 19901 ; iLesieur 199oh in MHD 
turbulence. Our simulation results described in Fig. (1) are 
fully consistent with this scenario. 



5 ANISOTROPIC TURBULENT CASCADES 

We find that the presence of a large scale magnetic field 
Bo (along the ^-direction) introduces anisotropy in the dis- 
tribution of energy in wavevector space such that the rms 
wavenumbers along (kit) and across (k±) the mean magnetic 
field Bo show a discrepancy i.e. kit ^ k± [see Fig. (3)]. We 



Figure 3. (Left) Evolution of fcii/fcj_ in anisotropic magnctofluid 
turbulence shows that kj_ becomes progressively dominant. 
(Right) Anisotropic magnetic field spectra along and across the 
Bo are consistent with the left panel. The spectra are computed 
in the steady state (close to 11 — 13 lo/vo) and are from the 
same simulation. The forcing mode spans the wavenumber band 
3 < k f < 5. 



Figure 4. Spectral distribution of the ratio |_B(fcjJ| 2 /|B(fc||)| 2 



employ the following diagnostics to monitor the evolution 
of the rms wavenumbers k± and k\\ in time. The rms k± 
mode is determined by averaging over the entire turbulent 
spectrum weighted by k±, thus 



k ± = (k±(t)) = 



\k±B(k,t)\ 2 



E k \B(k,tW 



Here (• • •) represents an average over the entire Fourier 
spectrum, k±_ — yfcl + fcf- Similarly, the evolution of the 
fey = k z (along the Bo direction) mode is determined by the 
following relation, 
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fc|, = <*li(t)) 



£Jb(m)I 2 



The wavenumbers k x , k y and k z are respectively along the 
x, y and z directions. It is clear from these expressions that 
the k± and feii modes exhibit isotropy when k± ~ feii . Any de- 
viation from this equality corresponds to spectral anisotropy. 
We follow the evolution of k± and fcii in our simulations. We 
find a disparity in the magnetic field fluctuation spectrum 
along 

|B(fc|l)| 2 = ^|B(fcx,fc||)j 2 dfcx, 

k ± 

and across 

|B(fcx)| 2 =^|i?(fcx,fc||)| 2 dfc|| 

the mean magnetic field (see Fig. 3, right panel). Note that 
the discrete summations } lf _ and ^2 k± are carried over 
fcii and k± modes respectively. The presence of the mean 
magnetic field inhibits turbulent cascades in the parallel (to 
Bo) direction and hence the characteristic modes (ku) along 
the mean magnetic field are suppressed, while the modes 
in the orthogonal direction remain unaffected. The evolu- 
tion in Fig. (3) therefore shows that the initial isotropic 
ratio fen ~ 1 progressively evolves towards anisotropy 
fe||/Jfc± < 1. 

The suppression of the ku mode is caused by the excita- 
tion of Alfven waves, which act to weaken spectral transfer 
along the direction of propagation. This can be understood 
as follows; We assume that the spectral transfer, mediated 
by propagating Alfven waves, can be described by a three 
wave interaction mechanism, for which the frequency and 
wa venumber resonance criteria are, respective ly, expressed 
by (|Shebalin et al. 19831 ; iMatthaeus et al 19981 ) 

±^3 = U)\ — UJ'2, 

and 

k 3 = ki +k 2 . 

The frequency-wavenumber resonance conditions indicate 
that two Alfven waves (a>i,ki) and (u^jks) mutually inter- 
act and give rise to the third wave (a>3,ka). Such conditions 
could, in principle, hold for a set of infinite waves as the 
indices '1' and '2' are merely dummy indices. With the help 
of the Alfven wave dispersion relation (iv = k ■ Bo) and the 
component of wavenumber matching relation along Bo, i.e. 
fcari = feln + fca ii , we infer that either the ki,, = or fc2„ = 0. 
It follows from the frequency-wavenumber resonance con- 
ditions that either the ki,, = or fca = 0. Owing to the 
absence of one of the parallel (to Bo) components of the 
modes, nonlinear mode coupling interaction becomes ineffi- 
cient in transferring the inertial range spectral energy in the 
parallel direction. Hence there is very little cascading along 
the magnetic field direction. Thus, the parallel wavenum- 
bers (fey) appear to be suppressed and the spectral cascade 
mainly occurs in the perpendicular wavenumbers (k±). Con- 
sequently, the magnetic field spectrum along the mean Bo is 
depleted. This, we suggest, explains the wavenumber dispar- 
ity fey 7^ k± observed in our simulations [see Fig 3]. While 
the turbulent cascades are effected locally by the presence 



of the mean magnetic field, the 3D volume averaged inertial 
range spectra in our simulations continue to exhibit a power 
law close to that of Fig. ((TJ). 

The question of anisotropic inertial range cas- 
cades has long been debated in the context of MHD 
turbule nce. It dates ba ck to the seminal work of Irosh- 
nikov ijlroshnikov 1 963) and Kraichnan (1965), who 
first pointed out that the presence of a large-scale or 
self-consistently generated magnetic field influences the 
spec tral power cascade mechanism in a complicated man- 
ner l|lroshnikov 19631 ; iKraichnan 1 965). Kraichnan (1965) 
addressed the interaction of magnetized turbulent eddies 
with Alfven waves, excited as a result of a mean magnetic 
field. Owing to the presence of waves in a MHD fluid, 
turbulent correlations between velocity and magnetic field 
and the corresponding energy transfer time are determined 
primarily by r ~ (bok)~ , where bo is a typical amplitude 
of a local magnetic field and is dimensionally identical to 
that of the velocity field by virtue of Elsasser's symmetry 
(Kraichnan 1965). According to Kraichnan (1965), it is 
this time scale that leads to the modification in the energy 
transfer (or cascade) because of the wave-turbulent eddy 
interaction process and without this modification the 
energy cascade rates would be determined by typical hydro- 
dynamic eddy interaction time. The ass umption of isotropy 
in the Kolmogorov energy spectrum l|Kolmogorov 1941) 
was later modifie d to include the mean or large scale 
magnetic f ield (IKraichnan 19651; IShebalin et al. 1983 ; 
Higdon 19841; iHigdon 198rj ; iGoldreich fc Sridhar 19951 ; 
Matthaeus et al 19981 ) . in part, because the presence of 
waves in MHD turbulence can significantly influence 
the energy cascade dynamics in the wavenumber space. 
An attractive analy sis was presented by Shebalin et al 
(|Shebalin et al. 19831 ). who demonstrated that the presence 
of a mean (or dc) magnetic field introduces a spectral 
anisotropy in the MHD turbulent spectrum. The observed 
anisotropic cascade in their work was understood to be due 
to the presence of Alfven waves, which are excited and prop- 
agate along the mean magnetic field. The non-dispersive 
propagating Alfven waves in the presence of a mean mag- 
netic field give rise to distinct energy cascade rates along 
and across the large-scale magnetic field, thereby leading 
to an asymmetry in the spectral transfer rates. This can be 
understood as follows; The inertial range energy cascade 
rates depend on wavenumber (or modes). Since the presence 
of the background magnetic field depletes the parallel mode 
(i.e. fcii or k z ) and leaves the perpendicular mode k± 
unaffected, the spectral energy transfer corresponding to 
the fcii mode is suppressed in the direction of the mean or 
background magnetic field. By contrast, the perpendicular 
transfer of spectral energy remains unaffected. This is 
quantified by Fig (3), right panel in which the spectrum 
of |B(fc||)| 2 is suppressed as compared to that of |B(fex)| 2 - 
It is because of this disparity, the energy cascades along 
and across the background magnetic field are distinct. 
Such asymmetric cascades, in agreement with arguments 
based on the frequency-wavenumber resonance conditions 
described a s above, produce spectral anisotropy in M HD 
turbulence (|Shebalin et al. 19831 ; IMatthaeus et al 19981 ). In 
the context of the solar wind plasma, fast (> 500km/s) and 
slow (< 400km/s) streams are dominated respe ctively by 
fcii and k± anisotropic modes (|Dasso et al 20051 ). A more 
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quantitative treatment of anisot ropic cascade was put for - 
ward by Goldriech and Sridhar l|Goldreich fc Sridhar 199a ) 
for MHD turbulence, suggesting that the Kolmogorov 
energy spectrum corresponds to wavenumbers perpendic- 
ular to the local magnetic field so that E(k±) ~ fc^ J,/3 
(where k± is wavenumber perpendicular to the local 
magnetic field), whereas the parallel wave number scales 
as fell oc k 2 / 3 . Boldyrev (2005) proposes that a G oldreich- 
Sridhar like spectrum l|Goldreich fc Sridhar 19951 ) results 
from a weak large-scale magnetic field, while the limit 
of strong anisotropy, that is, strong large-scale mag- 
netic field, corresponds to th e froshnikov-Kraichnan 
ijlroshnikov 19631 ; iKraichnan 19 65) scaling of the spectrum. 
This suggests that there exists an asymmetry in the energy 
cascade mechanism, which is believed to be due primarily 
to the presence of large-scale local magnetic fields in the 
turbulent MHD flow. Our simulation results describing the 
anisotropic MHD turbulent cascades, depicted in Fig. (3), 
are further consistent with (jShebalin et al. 19831 ). 



6 SUMMARY 

Our work proposes a self-consistent physical paradigm 
for the development of the density fluctuations spectrum 
in solar wind plasma in the context of compressible, 
driven-dissipative, anisotropic 3D MHD turbulence. Under- 
standing turbulent cascades in the MHD plasma is critical 
to many astrophysical phenomena. These range from 
understanding the role of waves and nonlinear cascades in 
the evolution of the solar wind, structure formation at the 
largest scales, cosmic ray scattering and energization by 
solar wind turbulence at the smallest scales and the heating 
of the solar wind (Scalo & Elmegreen 2004, Elmegreen 
& Scalo 2004, Goldreich & Sridhar 1995, Zank 1999, 
Goldstein et al 1995, Scalo et al 1998, Goldreich 2001) to 



ISM jRov et al 20081; Havcrkorn 20081; Rvu et al 2008 


Rickett 20071: IWillett 2005: 


Elmeereen 199S 


Dickev et al 200 ll iMinter & Spanrier 1996 


)■ 



We find from our 3D (decaying turbulence) simula- 
tions that a fc _5//3 density fluctuation spectrum emerges in 
fully developed compressible MHD turbulence from nonlin- 
ear mode coupling interactions that lead to the migration of 
spectral energy in the transverse (i.e. k _L U) Alfvenic fluc- 
tuations, while the longitudinal "compressional modes" cor- 
responding to k || U fluctuations make an insignificant con- 
tribution to the spectral transfer of inertial range turbulent 
energy. The explanation, in part, resides with the evolution- 
ary characteristics of the MHD plasma that governs the evo- 
lution of the non-solenoidal velocity field in the momentum 
field. It is the non-solenoidal component of plasma motions 
that describes the high frequency contribution correspond- 
ing to the acoustic time-scales in the modified pseudosound 
relationship (Montgomery et al 1987; Matthaeus et al 1988; 
Zank & Matthaeus 1990; Zank & Matthaeus 1993). What 
is notable in our present work is we find a self-consistent 
evolution of a Kolmogorov-like density fluctuation spectrum 
in MHD turbulence that results primarily from turbulent 
damping of non-solenoidal modes that constitute fast and 
slow propagating magnetoacoustic compressional perturba- 
tions. These are essentially a higher frequency (compared 



with the Alfvenic waves) component that evolve on acoustic 
time-scales and can lead to a "pseudosound relationship" as 
identified in the nearly incompressible theory (Matthaeus et 
al 1988; Zank & Matthaeus 1990; Zank & Matthaeus 1993; 
Bayly et al; Shaikh & Zank 2004a,b,c, 2006, 2007). The most 
significant point to emerge from our simulation is the dimin- 
ishing of the high frequency component that is related to the 
damping of compressible plasma motion. This further leads 
to the dissipation of the small scale and high frequency com- 
pressive turbulent modes. Consequently, the MHD plasma 
relaxes towards a nearly incompressible state where the den- 
sity is convected passively by the velocity field and eventu- 
ally develops a fc -5 / 3 spectrum. This physical picture sug- 
gests that a nearly incompressible state develops naturally 
from a compressive MHD magnetoplasma in the solar wind. 

The higher resolution in-situ solar wind velocity field 
observations at 1AU show a typical spectral index in the 
range from 1.4 to 1.5 (Podesta et al 2007; Tessein et al. 
2009). The latter clearly differs from our simulations in 
which the spectrum of density field is close to 1.67. To 
this end, we should point out that the solar wind obser- 
vations at 1 AU contradict the fundamental theme of our 
model which states that the density fluctuations are con- 
vected passively by the velocity field and thereby acquire 
similar spectral power-law. Owing thus to the differences 
between the observations and our simulations, our model de- 
scribing the passive convection of the density field through 
the nearly incompressible velocity field may not hold near 
1 AU. Consequently, the solar wind density fluctuations 
(~ fc~ 5 / 3 ) at 1 AU cannot be correlated with the velocity 
field (~ fc _3//2 ) through the passive-advection phenomenol- 
ogy. We believe that the deviation of the density and ve- 
locity fields, leading thereby to the obvious differences near 
1 AU, may arise from a number of physical processes, as 
noted below. Firstly, the physical processes driving the ve- 
locity spectrum near 1 AU are different from those beyond 
1 AU. Therefore, the evolutionary characteristics associated 
with the velocity spectrum near 1 AU (describing a 3/2-like 
spectrum) can certainly not be representive of the distant 
outer heliospheric turbulence spectrum. Secondly, close to 
1 AU, stream interactions, shear instabilities, compressional 
modes, and many other physical processes can plausibly al- 
ter the velocity field spectrum. It is unclear from the work 
of Podesta et al (2007) and Tessien et al (2009) whether 
any of these processes have a significant influence on the ob- 
served spectral indices. To further clarify this point, we have 
carried out more simulations to distinguish the evolution of 
kinetic energy, and incompressible, compressible, and total 
velocity spectra (not shown in this paper). We find that 
the compressive velocity field exhibits a flatter spectrum, 
while the incompressible velocity field follows a 5/3 spec- 
trum. This could very well mean that the velocity field, if 
is, driven or dominated by compressive modes, has a spec- 
trum may be a flatter than 5/3. Thirdly, the Alfvenic inter- 
actions, often invoked to explain the descrepnacy between 
the Kolmogorov-like (5/3) and Kraichnan-like (3/2) spectra, 
are ascribed to Alfven waves in MHD turbulence (Irosh- 
nikov 1963; Kraichnan 1965; Shebalin et al 1983; Biskamp 
2003). The latter inhibits energy cascades along the direc- 
tion of propagation. The spectral transfer of inertial range 
turbulent energy is therefore suppressed along the mean 
Bo, whereas the perpendicular cascade is governed predomi- 
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nantly by the bydrodynamic-like interactions. Consequently, 
the energy spectrum is dominated by hydrodynamic-like 
processes which lead to a Kolmogorov-like (5/3) inertial 
range turbulent spectrum. Hence a Kolmogorov-like (5/3) 
spectrum emerges in MHD turbulence when the nonlinear 
interactions are dominated by hydrodynamic-like eddies. By 
contrast, magnetic field eddies governing the Alfvenic inter- 
actions lead to a Kraichnan-like (3/2) spectrum in MHD 
turbulence. Thus, the disparate time scales associated with 
Alfvenic and compressive modes may flatten out the velocity 
field but not the density field. Our model, thus relating the 
density to the velocity field spectrum, may not hold near 1 
AU in the circumstances where the velocity fluctuations are 
driven by the processes as described above. 

The Kolmogorov-like fc -5 ^ 3 spectrum in density, veloc- 
ity and magnetic fields resulting from our simulation (see Fig 
1) are fully consistent with those of Tilley & Pudritz (2007), 
Mac Low et al (1998, 1989), Biskamp (2003), Padoan & 
Nordlund (1999), Stone et al (1998), Goldstein et al (1995), 
Goldreich & Sridhar (1995) and with others that are de- 
scribed elsewhere in our paper. 
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